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In a recent experiment, Young et al. [Nature 505, 528 (2014)] observed a metal to insulator transition as 
well as transport through helical edge states in monolayer graphene under a strong, tilted magnetic field. Under 
such conditions, the bulk is a magnetic insulator which can exhibit metallic conduction through helical edges. 
It was found that the two-terminal conductance of the helical channels deviates from the expected quantized 
value (= e^//i per edge, at zero temperature). Motivated by this observation, we study the effect of disorder on 
the conduction through the edge channels. We show that, unlike for helical edges of topological insulators in 
semiconducting quantum wells, a disorder Rashba spin-orbit coupling does not lead to backscattering, at least 
to leading order. Instead, we find that the lack of perfect anti-alignment of the electron spins in the helical 
channels to be the most likely cause for backscattering arising from scalar (i.e. spin-independent) impurities. 
The intrinsic spin-orbit coupling and other time-reversal symmetry breaking and/or sublattice-parity breaking 
potentials also lead to (sub-leading) corrections to the channel conductance. 


I. INTRODUCTION 

The quantum spin hall effect (QSHEjHH is an interesting 
phenomenon related to the emergence of an insulating elec¬ 
tronic state in the bulk of a two-dimensional material whose 
one-dimensional edge remains metallic and exhibits quan¬ 
tized conductance at zero termerature. As first discussed by 
Kane and Mele for graphene,[3 in the QSHE state, the helical 
edges can transport a spin current without dissipation as long 
as time-reversal symmetry is not broken. Unfortunately, the 
spin-orbit coupling is very small in graphend®^ (~ 10 peV) 
and therefore the effect is not experimentally accessible in this 
material. 

Nevertheless, in a recent experiment!^ transport through 
helical edges has been observed in monolayer graphene when 
the latter is submitted to a strong, tilted magnetic field 
B = (i?||,0, B±). Under such conditions, the bulk be- 



Figure 1. Schematic representation of the system considered in this 
work: A large metallic cluster (black dot) located near an armchair 
graphene edge induces potential and spin-orbit coupling by proxim¬ 
ity. The size of the cluster is taken to be comparable to the magnetic 
length (Rm ~ Ib, not drawn to scale). The cluster induces backscat¬ 
tering between the two counter propagating helical edge channels. 


comes a magnetic insulator whilst the edge can remain metal¬ 
lic. By increasing the strength of the component of B par¬ 
allel to the graphene plane (i?||), the magnetic state can be 
tuned from anti-ferromagnet (AF), to canted anti-ferromagnet 
(CAE), an d fina lly to ferromagnet (FM).E2ffi31 As shown by 
Kharitonov,(see also Ref. [fSl) the nature of the magnetic 
insulating is determined by the competition of the Zeeman and 
interaction energies. At the edge of the magnetic insulator in 
the CAE and EM states, electrons propagate in opposite di¬ 
rections with ( roughly) opposite spins, that is, a helical edge 
is formed.Esniltni Although time-reversal symmetry is broken 
explicitly by the external magnetic field, the two-terminal con¬ 
ductance of the the helical edge approaches the quantized con¬ 
ductance at the largest magnetic field: The conductance per 
edge is G = Go ~ SG, where Go = jh and SG/Go 0.1 
at the lowest accessible temperatures in the EM state.1121 

The physical origin of deviation 5G is not entirely clear 
and it is the main purpose of this work to investigate how dis¬ 
order can contribute to it. To this end, we assume that a po¬ 
tential scatterer is located near an otherwise perfect armchair 
graphene edge (see Fig. [^. Our model is not intended to be 
a faithful description of the experimental situation found in 
Ref [121 but we believe that it serves well to the purpose of 
investigating the effects of disorder-induced dissipation at the 
helical edges. However, if necessary, the model can be eas¬ 
ily generalized to extended disorder potentials, a task that we 
shall not pursue here. 

Nevertheless, the situation envisaged here can exist in an 
actual experimental device, where the scatterer may be a con- 
tacP^or a metallic cluster located near the edge of the device, 
for instance. If the cluster/contact contains heavy atoms (e.g. 
golcP^, the latter can locally induce by proximity a spin-orbit 
coupling (SOC) as well as providing a source for potential 
scattering.E^Jl^ Therefore, one of the important goals of this 
work is to investigate the role of the SOC-induced dissipation 
at the edge. Indeed, in the case of semiconducting quantum 
wells, it is known that Rashba-like disorder SOC has been 
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identified as one of the possible causes for backs catteringP. 
In addition, it was suggested by the authors of Ref. [12] that 
Rashba-like disorder could be an explanation for the absence 
of perfect quantization of the edge conductance. 

In this article, we derive the low-energy effective descrip¬ 
tion of the above model accounting for electron-electron in¬ 
teractions in the helical edge channels, which are crucial in 
determining the temperature dependence of 6G. The effective 
low-energy model turns out to be a version of the Kane-Fisher 
model for an impurity in a Tomonaga-Luttinger liquid.l^ Our 
analysis shows that the leading contributions to 6G are po¬ 
tential scattering and, to a lesser extent, the so-called intrinsic 
(or Kane-Mele) type SOC. However, we hnd that, to leading 
order, the Rashba SOC does not lead to backscattering. In ad¬ 
dition, our analysis indicates that other sources of backscatter¬ 
ing sources are potentials that explicitly violate the sublattice 
inter-change (parity) symmetry and/or time-reversal symme¬ 
try. 

The rest of this article is c^anized as follows. In Sec. |I^ we 
review Kharitonov’s results^lfor the the energy dispersion of 
the armchair edge. Hence, we derive the effective low energy 
model for the edge states, which turns out to be, accounting 
for electron-electron interaction, a Tomonaga-Luttinger liquid 
with two counter-propagating edge modes of different Fermi 
velocity. In Sec. m we discuss the properties of a general 
scatterer, which may correspond to a large metal cluster lo¬ 
cated near the edge of the sample. In Sec. IV we obtain 
the finite temperature conductance corrections within linear 
response theory. A summary and a brief discussion of the 
possible extensions to this work are provided in Sec. |V] The 
most technical details of the calculation have been relegated 
to the Appendices. In what follows, we work in units where 
H = kg = i- 


II. LOW-ENERGY DESCRIPTION OE THE EDGE 

In order to make the article more self-contained, we will re¬ 
view the main results of Kharitonov’s calculation for a semi- 
inhnite graphene layer submitted to a strong, tilted magnetic 
field. Within the k ■ p continuum description, the low-energy 
effective Hamiltonian of graphene under a tilted magnetic 
field can be written as a sum of three terms, H = K+V+H^, 
where 

K = J dr'ii^r){vFTz(J:^nx+VF(yyny)'i/{r), (1) 

If 

V = - y drdr''k^(r)'I'^(r')-j-^T'(r')T'(r), (2) 

Hz = -]^gsk-B J dr T'^(r)(s • B)T'(r); (3) 



Figure 2. Low energy band structure for non-interacting electrons in 
graphene terminated at an armchair edge as af function of pi b where 
p is the momentum along the edge (the edge is located at pis ~ 0) 
and IB the magnetic length. Electrons occupy the region pis ^ 0. 
The energy is in \/2 vp jlB units. The zero Landau level (E — 0 for 
pIb —> -foo) consists of one hole and one electron band, which are 
split near the edge. 


(e < 0 is the electron charge). The matrices (Tq, (tq.) act 
upon the sublattice (valley) pseudo-spin and Sq act upon the 
spin degree of freedom. Note that the z-axis is perpendicu¬ 
lar to the graphene plane and the x-axis will be taken along 
the armchair edge). The Zeeman energy is proportional to the 
total magnetic field B — (if||, 0, B±), where B± (i3||) is the 
magnetic field along the perpendicular (parallel) direction to 
graphene plane. 

Since the largest energy scale of quantum hall system is set 
by the cyclotron energy, ojc = vf/Ib where Ib is the mag¬ 
netic length, it is useful to neglect the Zeeman and interac¬ 
tion terms in the Hamiltonian and diagonalize the kinetic term 
first. Furthermore, in the study of edge effects, two kinds of 
edge terminations are commonly used, namely the armchair 
and zig-zag edge.l^SEl in this work, for the sake of simplic¬ 
ity, we use armchair edge, which we assume to be located at 
y = Q. The energy spectrum of the kinetic energy part of 
the Hamiltonian for a semi-infinite armchair edge is shown in 
Fig. 0 Note that, within this approximation, all the energy 
levels spin-degenerate and therefore the ground state is non¬ 
magnetic. 

For what follows, it is useful to express the fermionic field 
operators 'k (r) and (r) in the basis of Landau level orbitals 
that diagonalize the kinetic energy K: 


K is the kinetic energy, V the Coulomb interaction, and 
Hz the Zeeman term; the couplings ps, gs and B are the 
Lande g-factor, the Bohr magneton and the magnitude of the 
applied (total) magnetic field, respectively. The Fermi veloc¬ 
ity of electrons in graphene is vf ~ c/300 (c being the speed 
of light), and n = p — eA{r)lc is the kinetic momentum 


'^{r) = (j)lp^{r)Cnprs- (4) 

nprs 

where p is the momentum along the edge (i.e. x-axis) and 
n = 0, ±1, ±2,... is the Landau level quantum number. In 
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the Landau gauge where A = {B±y, 0,0), 



^zK'r 
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0,p) 
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(5) 

( 6 ) 

(7) 

( 8 ) 


for n = 0. The states \n,p) are the solutions of the Lan¬ 
dau level problem of the Schrodinger equation in the Landau 
gauge. In particular, 


(r|n = 0,p) 


gipx exp 21 ^^ j 


(9) 


where Ib = ^ is the magnetic length and the length 
of the (armchair) edge. Finally, the spinors: 


»?T = 




( 10 ) 


Note that in the n = 0 Landau level (OLL), the sublat¬ 
tice and valley (pseudo-)spins are aligned. This feature will 
probe extremely important in the discussion that follows. Ac¬ 
counting for the spin degeneracy, each single-particle state in 
the OLL is four-fold degenerate corresponding to s =t, i and 
T = K, K'. 


A. Zeeman and interaction energy effects 


In order to fully determine the ground state, we next take 
into account the Coulomb interaction and the Zeeman term. 
In neutral (i.e. undoped) graphene, the Landau level of zero 
energy (OLL) is half filled and therefore all Landau levels be¬ 
low (above) OLL are completely filled (empty). At tempera¬ 
tures far below the cyclotron energy (i.e. T ^ Wc) and close 
to half-filling, the low energy dynamics is determined by the 
electrons in the OLL. As a result, we project the Coulomb 
interaction (Eq. onto the subspace of states belonging to 
the OLL. Thus, in what follows, we shall suppress the Landau 
level index n as it is implicitly understood that n = 0. In addi¬ 
tion, we will suppress the sublattice index because in the OLL 
the pseudo-spin and valley spin are aligned. The Coulomb in¬ 
teraction in the OLL subspace is therefore formally described 
bjEH 


VoLL = ^ H 

ot—0,x^y,z pi...p4 

'-piO's('^a)a'o-'Cp2cr's j ( 11 ) 

s (T,cr' 


where : ... : stands for normal ordering. The Coulomb am¬ 
plitudes Ua have been discussed in Refs. [13] and IZTl and we 
will not dwell on them here. Nevertheless, it is worth men¬ 
tioning that the Coulomb amplitudes that do not conserve val¬ 
ley quantum number are exponentially suppressed by an ex¬ 
ponential factor ~ Yvbere a is the distance between 

the carbon atoms in graphene. Thus, the intra-valley ampli¬ 
tude, uq, is dominant one.^^ In fact, if only uq is retained in 
Eq. ([TTJ, the Hamiltonian and the ground state would exhibit 
SU(4) symmetry (see ReflZTland references therein). 

However, when the sub-leading amplitudes Mq=i, 2,3 that 
violate the SU(4) symmetry are also included, together with 
the Zeeman energ y, a qua ntum hall magnet can be stabilized 
as the ground stateP^EHIlits magnetic order can be uncovered 
by relying upon mean-field the ory f or which the order param¬ 
eter is This order parameter can 

be written as follows: 

P = Xaxi+Xbxl, (12) 

where Xa = \K) <8 |S'a) and Xb = \K') 0 \Sb) are the trial 
spin-valley spinors. For the semi-infinite system, we shall ne¬ 
glect any deviatio n of the order parameter P from its bulk 
value near the edge.EE® We shall also neglect the dynamics 
of the order parameter, which can lead to additional dissipa¬ 
tion mechanisms in clean systems (via coupling with the bulk 
Goldstone modes), different from those discussed here.li^ In 
spin space, the spinor accepts the following parametrization: 


\Sa,b) 


f cos ^ 

y±sin 



(13) 


where Og is the angle between the spin polarization and the 
total magnetic field B-, (pg is the azimuthal angle around the 
total magnetic field . The positive sign applies to Sa and the 
negative sign to 5*^1 13 1 141 above parametrization assumes 
that the spin-quantization axis points along the direction of the 
total magnetic field B. Note that this choice is different from 
the one assumed for the operator s^, which measures the spin- 
projection perpendicular to the plane of graphene and Sx, Sy, 
which measure the in-plane spin projection. Both sets of spin 
operators Sq, and (a = x, y, z) are related by a rotation: 


= Sz COS0-I-Sa;Sin0, (14) 

sinS-f Sj, COS0, (15) 

Sy = Sy. (16) 


By minimizing the kinetic energy and Zeeman term to¬ 
gether with the Coulomb interaction within the mean-field ap¬ 
proach, the single-particle spectrum shown in Fig. is ob- 
tained.Ii^The Zeeman term in the Hamiltonian lifts the degen¬ 
eracy between the spin up and spin down bands. However, 
due to the exchange part of the Coulomb interaction, the OLL 
bands near the edge anti-cross (for pis ~ 0.25) and develop 
an energy gap. The magnitude of the gap is given by,G^ 



(CAF) 

(F) 


(17) 
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Figure 3. Lowest energy band structure of a graphene armchair edge within mean field theory for the canted anti-ferromagnetic (CAF, left 
panel) and ferromagnetic (FM, right panel), respectively. The edge is located at pis « 0. As a result of the Zeeman and the Coulomb 
interaction, an energy gap develops in the hand structure, which is minimum near the edge. In electron doped graphene, the low energy 
dynamics are governed by the degrees of freedom near the Fermi points where the upper band intersects with the Fermi energy Ef (for hole 
doping Ef is within the lower band). 


where 6 = tan B^) is angle subtended by B with the 

z-axis perpendicular to the graphene plane, ez = \gs^^BB 
is the Zeeman energy, and A = uq -F 2|ui| + M 3 ; 61cr = 
cos“^ (ez/ 2 |M_L|) is the critical angle for which a quantum 
phase transition from the canted anti-ferromagnet (CAF) and 
ferromagnet (FM) takes place when the Zeeman energy takes 
over the exchange energy. 


In the CAF state, when the chemical potential lies above 
(below) the gap in the upper (lower) band, two metallic chan¬ 
nels appear at the edge (corresponding to the two Fermi points 
in Fig.j^where the energy band intersects the chemical poten¬ 
tial). The corresponding states are not eigenstates of = 
s ■ B/B and are not mutually orthogonal. This has a signifi¬ 
cant effect on the suppression of conductance from its perfect 


quantized value Go, as we will see in Sec. Ill 


As the angle 9 is tuned across its critical value, 0cr, and the 
system becomes a FM, the two conducting channels along the 
edge still exists. Unlike in the CAF, they are eigenstates of 
(cf. Eq. [I4| i and therefore mutually orthogonal. Nevertheless, 
we note that this is the result of the approximation made in 
Ref [ 13 ] which neglects the variation of the order parameter 
near the edge. The more sophisticated treatment provided in 
Ref [ 15 ] indicates that spin orientation of the edge states re¬ 
mains slight canted in the FM, only becoming perfectly anti¬ 
aligned at large values of i?jj. However, to the extent they 
can be assumed to be perfectly spin anti-aligned, the FM edge 
channels can transport spin and charge current without dis¬ 
sipation if there were no impurities. However, different from 
the QSHE in semiconducting quantum wells, the QSHE in the 
EM quantum hall insulator is not protected by time-reversal 
symmetry but the U(l) symmetry generated by the total s^. 


energy Hamiltonian for the edge modes can be obtained; 

■ ■, (18) 

c=R,L p 


where 'ipR^p) ~ Cp^ ^^+p and V’l(p) = Cp^L-p for \p\ < A 
and A ^ Ep/ max{vii, vp} and 


vr = 


dE{p) 


dp 


VL = 


dE{p) 


dp 


(19) 


The creation and annihilation operator obey the usual fermion 
anti-commutation relation {'0c(p), V'c'( fOI = Spy6c,p 
where c = {R,L). Note that, unlike the QSHE in semicon¬ 
ductor quantum wells for which the states at the two channels 
are a Kramer’s pair,l^ in the present case the right and left 
movers can have different Eermi velocities (i.e. vr ^ vp). 

In addition, since the degrees of freedom near the edge are 
gapless, we need to take into account the effect of electron- 
electron interactions. In experimental conditions, due to the 
presence of a nearby metallic gate used to tune the Eermi 
energy Ep, the Coulomb interaction between electrons is 
screened to yield a short-range interaction: 


lAdge — ‘^'^ 9 RL 


dx : 'ilj^p^{x)ipR[x)'tl;\{x)'tl;p{x) : . (20) 


Note that I4dge is not necessarily repulsive as the Coulomb 
interaction may be over screened at low energies and become 
an effective attractive interaction. As a consequence, below 
we have considered both signs for gpp. Thus, the complete 
low-energy Hamiltonian for a clean armchair edge reads: 

Hedge = Hq + t4dge- (21) 


III. SCATTERING POTENTIALS 


In both the CAE and EM states, at T ^ minjAedge, Ep}, 
the low-energy physics is controlled by the fermionic degrees 
of freedom near the Eermi energy Ep. Upon linearlizing the 
band around the Eermi points defined by the intersects {pp,p 
and pp,r) of the bands with Ep (see Eig. [^, an effective low- 


A. General form of a localized scatterer potential 

We next analyse how the presence of a scatterer affects con¬ 
duction along the edge. To this end, it is instructive to consider 
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the most general form of the potential created by a generic 
scatterer located near the edge of the sample. Generally, the 
latter can be written as follows; 

Vm = j dr ^'^(r)y(r)r['(r), (22) 

where V (r) is a 8 x 8 potential matrix acting on both sub¬ 
lattice, valley, and spin degrees of freedom. Since a tilted 
magnetic field break most symmetries of the graphene Hamil¬ 
tonian, the most general form of V (r) is a linear combination 
of the 64 matrices TaCrps^, where a, /3 ,7 = 0, a;, y, z (0 refer¬ 
ring to the identity matrix): 

V{r) = Y, (23) 

ck,/3,7 

where the relative strength of each potential Vai 3 .y{r) depends 
on the microscopic details of the scatterer, the strength of the 
external magnetic field B± and By, as well as other pertur¬ 
bations (boundary conditions at the edge). Therefore, Vm 
will contain terms that break the symmetries of the zero-field 
Hamiltonian. Those terms are generated by the projection 
onto the low-energy subspace containing the zero energy Lan¬ 
dau levels. In what follows, we shall discuss which terms are 
expected to yield a dominant contribution to Vm- Rather than 
trying to be exhaustive, our choice will be physically moti¬ 
vated. However, we shall first take a detour to discuss some 
general properties of the OLL projection. 


B. Projection of sublattice-valley Operators onto the OLL 

In order to obtain a low-energy effective description of the 
edge channels in the presence of the scatterer, its potential 
must be projected onto the OLL. To leading order, the scatterer 
potential in the OLL is given by 

VoLL = VV{r)V + ■■■ (24) 

where V is the projection operator in the OLL. The terms con¬ 
tained in the ellipsis represent corrections arising from virtual 
transitions to Landau levels with n ^ Q. Since the scatterer 
potential introduced above acts upon the sub-lattice, valley 
and spin degrees of freedom, it is worth dwelling a bit on 
the explicit form of the OLL projection operator V and how 
it affects the Pauli matrices TaCrps-y. 

The operator V can can be written as follows: 

V = Vsv®'Po, (25) 

Ps^ = \AK){AK\ + \BK'){BK'l (26) 

Vo=\n = 0,p){n = Q,p\ (27) 

The first term in the right hand-side of Eq. ( |25] l {Vsv) acts 
on the sub-lattice and valley pseudo-spin degrees of freedom, 
whereas the second term iVo) acts on the other orbital degrees 
of freedom. Since sub-lattice spin and valley spin are aligned 


in OLL, when projecting the Pauli matrices using Pav, we ob¬ 


tain the following: 

VsyT^a^Pav = VsvTz'J^'Psv = (28) 

'PsvT~a'^'Psv = , 'PsvT'''(T~'Psv = , (29) 

PsvTzCTzPsv = Vsv, PsvT^OqVsv = 0 , (30) 

VsvToa^Vsv = 0, VsvTzO^Vsv = 0, (31) 

PsvT^UzPsv = 0, 'PsvT^Cr'^Psv = 0, (32) 

PsvT~a~Vav = 0. (33) 

where we have introduced the operators; 

J:z = \AK){AK\-\BK'){BK'\, (34) 

E+ = \BK){AK'\, E" = \AK'){BK\. (35) 


Thus, {7^sj,,E^,E+,S } behave effectively as Pauli matrices 
in the subspace of OLL states. 


C. Time-Reversal Invariant Scattering Potentials 

Although the most general impurity potential should not re¬ 
spect time-reversal symmetry as the latter is broken by the 
external magnetic field, it is worth beginning our analysis 
by considering various time-reversal symmetric (TRS) po¬ 
tentials. The reason is that TRS potentials are present in 
the microscopic Hamiltonian of graphene with impurities at 
zero field. However, non-TRS potentials need to be gener¬ 
ated by means of virtual transitions to higher energy states in 
the presence of the time-reversal symmetry breaking pertur¬ 
bations like the magnetic field. Therefore, according to the 
discussion in previous section, TRS potentials are expected to 
have larger strength relative to the non-TRS ones. In addition, 
the discussion that follows is also motivated by the sugges¬ 
tion made in Ref. [T2| attributing the absence of conductance 
quantization to the existence of s catter ers proximity-induced 
Rashba spin-orbit coupling (SOC).l®3 

TRS potentials can be divided into two classes depending 
on how they act on the spin degree of freedom: Those propor¬ 
tional to which are termed ‘scalar’ potentials, and those 
proportional to which correspond to SOC. These two 

classes can be further subdivided into those inducing only 
intra-valley scattering and those inducing intervalley scatter¬ 
ing depending on whether they commute with or not. Note 
that, typically, the strength of the scalar potentials is much 
larger than the strength of SOC potential^^Si, 

An example for a TRS scalar potential is; 

Vs{r) = Vb(»’)so + Vx{r)Ta;So 

+Vxx{r)az:TxSo -I- (36) 

where the second and third terms [c>c 14 (r), 14a; (f)] and oth¬ 
ers similar to them not explicitly written describe interval¬ 
ley scattering. Strong inter-valley scattering requires that a 
rapid spatial variation of 14 (r) and 14a; (r) on the scale of 
the inter-carbon distance a in graphene. In other words, it 
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Scattering Operator 

OLL Projection 

ToCToSO 

sv So 

O'OTxSo 

0 

(^xTxSq 

So 

dzTzSz 

sv Sz 

^O'^z^z 

^2 S2: 

f^xTzSy 

0 

^y‘T~ 0 ^x 

0 


Table I. Projection of several types of time-reversal symmetirc scat¬ 
tering potentials onto the zero energy Landau Level of graphene 
(OLL). We have introduced the operator Sa, = (E"*" -|-E“)/2 


requires atomic scale disorder (i.e. point defects, edge rough¬ 
ness, etcetera). It is possible to estimate, within the Born ap¬ 
proximation, the reflection coefficient at the edge from such 
small-size scatterers, El which is ~ {e^/vF){a/lB)‘^■ How¬ 
ever, for the experimentally accessible values of B±, Is ~ 
26/y/B±{T) nm ^ a = 0.24 nm, implying that the interval¬ 
ley scattering terms in the above equation can be neglected. 
For the same reason, we shall also neglect the inter-valley 
SOC scattering terms below. 

Next, we consider the (intra-valley) SOC scattering poten¬ 
tials, which can take the following general form; 

14o(r) = Visoir)a^T;^Sz + Cpbiso(^’)r^Sz 

+ V-[ix{r)aa:TzSy -VRy{r)aySx, (37) 

where the first term on the right hand-side corresponds to 
the intrinsic (or Kane-Mele) type of SOC, and the third and 
fourth terms correspond to the Rashba-like SOC. Note that we 
have considered a generalized anisotropic Rashba, for which 
VRy(r). Such SOC potentials can be induced 
when the hexagonal point-group symmetry of graphene is 
broken down to a smaller symmetry group by proximity to 
substrates and clusters of heavy -metals.EE! The second term 
cx fpbiso(f) is absent in the microscopic Hamiltonian because 
it is forbidden by the parity symmetry that exchanges of the 
A and B sub-lattices. However, it can be generated by the ap¬ 
plication of parity-breaking fields like a tilted magnetic field 
or by the edge potential. Since this term is not present at the 
zero-field (for which parity symmetry is not broken), in a first 
approximation, we shall neglect it compared to the scalar and 
standard SOC terms. However, in the following section we 
shall see that this term can yield backscattering in the FM 
state. 

Upon projecting ( |36| ) and ( [J7| ) onto the OLL, we obtain; 
VVm^^V = f dr ^l{r) [Vo(r)so + ^so(t-)s2] 

(38) 

where 4'o(r) = J^prs (t^lpsi.r)cprs C^lir) = [4'o(r')]^) is the 
annihilation (creation) field operator in the OLL, respectively. 
Note that, in the above expression, measures the projec¬ 
tion of the electron spin along the axis perpendicular to the 


graphene layer and must be expressed in terms of and s/ 
(cf Eq. (fUji). 

It is noteworthy that Rashba SOC of the form given in 
Eq. @ gives no contribution to the OLL projected Vm 
(cf table [I]l. Therefore, it alone cannot suppress the con¬ 
ductance from its perfect quantized value, unlike the ca se of 
the QSHE in two-dimensional TRS topological insulators .EIH 
This is a direct consequence of the structure of OLL orbitals 
in graphene, for which the sub-lattice and valley pseudo-spins 
are locked. Landau level mixing arising from interactions, etc, 
can modify this conclusion slightly, but the corrections are ex¬ 
pected to be small and suppressed at the high magnetic fields. 
On the other hand, the scalar and intrinsic SOC potential are 
not affected by the OLL projection and therefore can lead to 
backscattering. 

To make further progress towards obtaining the effective 
Hamiltonian of a TRS scatterer at the edge, let us assume the 
size of the scatterer (e.g. a metal cluster of radius Rm) to 
be comparable to the magnetic length, Rm ^ Ib- Thus, we 
can approximate Vb(T') and UsoCt") by Dirac-delta functions, 
which, upon further projection on the low-energy modes near 
the edge, leads to 


Vm = qRm V'fl( 0 )i/'L( 0 ) + '0j,(O)'!/'fl(O) 


(39) 


where g has energy units. When combined with Eq. the 
Hamiltonian, H = Ffedge + Vm, is a version of the model 
studied by Kane and Eishei^lfor an impurity in a Tomonaga- 
Luttinger liquid. The only difference is that, in our model, 
right and left movers have different Eermi velocities. Note 
that, after projecting Vm on the edge, a term of the form 
b{x) has bee n drop ped since it does 
not contribute to the channel resistance The scattering po¬ 
tential coupling g contains contributions from both scalar and 
the intrinsic SOC potentials, i.e. g = gs + Qiso, where 


\/N{pr)N[pl) 


A{jpr)A{pl) , V 

^AF 


y/^V\so cos 9 
y/N{pR)N{pL) 


A[pr)A{pl) 

^AF 


(40) 

(41) 


with 


A{p) = Ap + ez - e{p) 

— \J^\f + [^(f) ~ + ez)]^ ; (42) 

e(p) is the single particle energy, which is displayed in Eig. 
and N{p) = y/A{pY + A\p is the normalization of the sin¬ 
gle particle wave function; 9 = tan“^ (FFy/iFj^) is the polar 
angle of the applied magnetic field, Ap = ^{uq — 2 |rti| + 
U 3 ) COS 9s and Aap = ^ (mq + 21 ui | -I- U 3 ) sin are the com¬ 
bination of the mean-field parameters that favour ferromag¬ 
netic and anti-ferromagnetic order, respectively. The relative 
strength of these potentials gs/Vo and giso/Vso has been plot¬ 
ted i n Fig . 1^ Note that gs/Vo 3> < 7 iso/biso- In addition, typi- 

caiiyli 2|201 

^ I^so> which means that the main contribution 
stems from the scalar part of the scatterer potential. 
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4. Normalized strength of the leading time-reversal invariant 
contributions to scatterer potential. The intrinsic SOC (dashed line), 
(?iso/Viso, and scalar potential (solid line), Qs/Vo, can be induced by 
proximity to a heavy metal cluster. Their strength is plotted against 
the angular deviation S6 — 9 — (in degrees); 0cr is the critical 
polar angle of the magnetic field for which quantum phase transition 
to ferromagnet (FM) occurs. The vanishing of the scattering potential 
strength in Ferromagnet (FM) is an artefact of the assumption that 
the scattering potential is entirely due to a combination of scalar and 
intrinsic (i.e. Kane-Mele type) spin-orbit coupling (SOC) along with 
the assumption of that the order parameter does not deviate from its 
bulk value. Note that, in the canted anti-ferromagnet corresponding 
to i9 < 0, the contribution of the scalar potential is larger than the 
SOC contribution. 
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D. Backscattering in the Ferromagnet 


in the FM. The reason is that, near the edge, the order param¬ 
eter can deviate from its bulk value and, as a consequence, the 
spins near the edge of the FM can be slightly canted.li^ This 
implies that the backscattering due to the scalar and intrinsic 
SOC potentials is suppressed relative to its values in the CAF 
but not they are completely eliminated. On the other hand, a 
Rashba SOC potential yields no backscattering to leading or¬ 
der. The latter result is independent of detailed form of the 
single-particle orbitals as it is entirely a consequence of the 
OLL projection. 

Nevertheless, it is also interesting to consider additional 
scattering potentials that have been neglected so far, which 
can also lead to backscattering. If we restrict ourselves to TRS 
potentials, as we did in the previous section, one possible can¬ 
didate, which by no means exhausts all the possibilities, is the 
following parity-breaking intrinsic SOC; 

^ J Vpbiso(f)^''''(T’)T2S2^'(r) H- (45) 

= J dr Fpbiso(f)^'(r)rz (s^ cos6» 

—sin0)'l'(r) + • • • (46) 

Other possible terms (included in the ellipsis) are proportional 
to TzSx,y Upon projection on the OLL, Eq. ( |46l l yields; 

VVZ^^V = J dr (4 cose 

sind) T'o(r) H- (47) 


As shown in Fig. the effective potential strength arising 
from the scalar potential and intrinsic SOC vanishes in the FM 
state. However, to understand the vanishing of gs,giso better, 
let us first recall the valley and spin structure of the single¬ 
particle orbitals in the FM state; 


\R)=\ 

l^\AK) + \BK')'^ 

1 ®lt), 

(43) 

\L)=\ 

flAK) - \BK')^ 

1 0 U). 

(44) 


Note that the second term on the right generates backscatter¬ 
ing in the FM whereas the first one, which commutes with 
does not. In the CAF, both terms lead to backscattering be¬ 
cause the spins are canted. However, as pointed out above, 
this term is not a standard type of SOC as it requires sub¬ 
lattice exchange (parity) symmetry to be broken. This may 
be caused by the substrate (BN), the tilted magnetic field, or 
the boundary conditions at the edge. However, estimating the 
strength of this term is beyond the scope of this work as this 
will require a detailed microscopic model of the system. 


Note that the in the FM the orbitals are eigenvectors of 
and therefore orthogonal. This means that no scalar (intra- or 
inter-valley) is effective in causing backscattering. Thus, only 
spin-dependent potentials can cause backscattering. However, 
the orbital part of the intrinsic SOC (^ ct^TzSz), cannot couple 
two edge states even if the spin operator does (recall that they 
are eigenstates of s^, not Sz)- 

However, we need to stress that the above conclusions rely 
on the orbital structure of the single-particle orbitals in the FM 
described by Eqs. ( |44l i and ( |4^. T he structure is the result of 
the assumption (made in RefTllA] and which we follow here) 
that the magnitude of the order parameter does not deviate 
from its bulk value. This approximation is expected to be ac- 
curatJi^for large values of Hy. However, a more sophisticated 
mean-field treatment like the one provided in Ref. [15] which 
allows for the spatial variation of the order parameter, would 
not automatically imply a vanishing impurity (back)scattering 


E. Non-time reversal symmetric scatterers 

Other possible potentials causing backscattering can break 
TRS. In particular, a potential of the form; 

= Y,dc.{r)s^, (48) 

Ct 

which describes the exchange interaction with a locally mag¬ 
netized scatterer (like a magnetized edgd^, gives a non-zero 
contribution to the backscattering potential the spins of the 
edge electrons are canted and therefore do not need to be com¬ 
pletely flipped in order to backscatter. 

However, at large values of i?||, for which canting is neg¬ 
ligible, will not produce backscattering in the EM 

state. The reason is that backscattering requires the potential 
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to act both on the sublattice-valley and spin degree of freedom. 
In other words, the potential must flip between (approximate) 
eigenstates of and meaning that the orbital part must be 
proportional and 'Ey = (E+ — E~)/2i. From the results 
in subsection IIIB we see that Ez arises from potentials pro¬ 
portional to aoTz or (T^To. Note that is odd under TRS, and 
since flipping the spin requires the potential be proportional 
to Sa, a term of the form TzSa is TRS. The latter results in 
the parity-symmetry breaking SOC discussed above. On the 
other hand, cr^ is TR even, and including the spin-part, leads 
to non-TRS potential of the form; 


We first consider the non-interacting case for which guL = 
0. The previous expression can be evaluated to yield: 

/ OO 7 

-G>{eV + u)G<{uj)\ (53) 

The non-interacting lesser and greater Green’s functions are 
G<{uj) = :^nF{uj) and G>(w) = ^(1 - npiu})), where 
np{uj) is the Fermi-Dirac distribution and c = R,L. Hence, 


j^NTRS 

*^2M 


KtrsMcTzSc 

a.—x,y,z 


(49) 


Since [V 2 ^^'®,Tz] = 0, this potential does not produce in¬ 
tervalley scattering. In addition, it also breaks sublattice ex¬ 
change parity symmetry like from Eq. ( |46| ). 

Finally, potentials containing Ey can arise from terms of the 
from CTqT /3 with a, /? = x,y, which describe inter-valley scat¬ 
tering and are therefore are suppressed by a factor ^ {a/lsY' 
as discussed above. 


IV. CORRECTIONS TO THE EDGE CONDUCTANCE 


In this section, we calculate the linear conductance of the 
edge channel in the presence of the scatterer. We focus on the 
doped CAF, for which the leading sources of backscattering 


have been identified in subsection III C and their strength can 


be estimated from their values in the zero-field Hamiltonian. 
The strength of the terms induced by TRS- and/or sublattice 
parity symmetry-breaking terms is much more difficult to esti¬ 
mate because they are generated upon integrating high-energy 
states and will also depend on the microscopic details of the 
scatterer. In what follows, we shall use perturbation theory to 
lowest order in the scattering potential strength, g, assuming 
that the latter is weak, i.e. g/oJc 1- 

As explained in Appendix |^the correction to the channel 
conductance within linear response theory is given by; 


5G = 


dls 

~dV 


v=o 


(50) 


where is the backscattered current. 


/ dNp \ 


(51) 


Np being the total number of right moving electrons. Within 
linear response theory, the (steady-state) backscattered current 


J361 


Ib =(2e) Im j C^At(i). 


(52) 


where the retarded correlation function is (f) = 

-i9{t){[A{t),A'^{0)]) and A = gRMi’piO)'ipLiO)- 


SG = 


(gRM)^ ^ 

-Go 

VrVl 


(54) 


For vp = vp, we recover the linear conductance of non¬ 
interacting ID channel with an impurity.l^ 

However, note that the conductance of the non-interacting 
edge is independent of the temperature, in contradiction with 
the the experimental observations!^ This problem can be 
solved by accounting for interaction effects (i.e. gpp ^ 0), 
which yield a temperature dependence of dG, in qualitative 
agreement with the experimental observation.^ Accounting 
for a non zero gpp is also necessary because interactions in 
ID have a drama tic effect on the stability of fermionic quasi¬ 
particles ) 34 | 37 | 38 | yyg calculate the correction to the con¬ 
ductance using the bosonization method, which allows us to 
treat interactions non-perturbatively (see Appendices). The 
result of this calculation reads: 


X B(iT,l-2i^)[V;(l-iF)-V'(^)]Go, (55) 


where K is the Luttinger parameter, which characterizes the 
strength of the interactions: K = \ for gpp = 0 and K < 1 
{K > 1) for repulsive (attractive) interactions. The parame¬ 
ters v± are the velocity of the eigenmodes (see Appendix 
for the definitions of K and u±); B{x,y) is the beta func¬ 
tion, V'(a^) is the digamma function. Note that the temperature 
dependence is a power-law, with an interaction dependent ex¬ 
ponent,!^ which vanishes in the non interacting case (i.e. for 
K = 1), in agreement with Eq. ( |54| ). 

In Fig.|^ we have plotted the two-terminal conductance per 
edge against the temperature (in cyclotron frequency units) for 
different values of Luttinger parameter K. Eor K > 1 (i.e. 
attractive interactions), the edge conductance approaches the 
quantized value as T —0. On the other hand, for K < 1 (re¬ 
pulsive interactions), the deviation from perfect quantization 
increases with decreasing T. For 5G ~ Gq linear response 
theory breaks down and we must rely on other methods, but 
this problem will not be analysed here. 

In Fig. 1^ the two-terminal conductance per edge is plotted 
against the deviation of the polar angle of the applied magnetic 
from the critical polar angle 9cr at constant temperature, for 
different Luttinger parameter K. 

Einally, let us briefly discuss the angular and temperature 
dependence of some of the parity-breaking SOC term dis¬ 
eased in Sec. Ill D This term is expected to yield a sizeable 
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Figure 5. (Color online) Two-terminal conductance per edge G = 
Go — SG (Go = fh is the quantized conductance) as a function 
of the temperature (in cyclotron frequency units), for different values 
of the Luttinger parameter K. This is plotted with the following 
parameters , = 45T, Bo_ = 5T, Rm = Ibs = 12nm. 
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Figure 6. (Color online) Two-terminal conductance per edge 
G = Go — SG (Go = I his the quantized conductance) as a func¬ 
tion of the polar angle of the applied magnetic field 59 = 6 — 9cr 
(the perpendicular magnetic field B±_ is fixed), for different values 
of the Luttinger parameter K. The angle 9cr = cos“^ [ez/(2|ui|)] 
is the critical angle for which the edge gap closes. Other calculation 
parameters used in this plot are chosen as follows T = 5 K, = 5 
T, Rm = Ib = 11.4 nm. 


contribution at large values of By in the FM state, when the 
cantin^Slof the edge channel spins is small. Upon projection 
onto the OLL, Eq. ( |46] l yields a term c>c sin 6, where 0 is the 
magnetic held polar angle. Hence, we obtain the following 
contribution of the conductance: 

<5Gpbiso(T, e) oc sin^ 9 (56) 

where T is the absolute temperature and K the Luttinger pa¬ 
rameter. However, since this term is not present in the zero- 
held Hamiltonian, we cannot estimate its strength, although 
we expected it to be small. 


V. SUMMARY AND OUTLOOK 

Using a simple single-scatterer model, we have discussed 
the effects of potential and spin-orbit disorder on conduc¬ 
tion through helical edges recently observed in monolayer 
graphene under a strong, tilted magnetic held.nS We have 
shown that the deviation from perfect quantization in the 
canted anti-ferromagnetic and ferromagnetic states can be 
understood in terms of potential (i.e. scalar) disorder and 
to a lesser extent spin-orbit disorder of the intrinsic type. 
Other types of disorder, such like sub-lattice exchange (par¬ 
ity) breaking spin-orbit coupling and time-reversal symmetry¬ 
breaking spin-orbit coupling terms are expected to yield 
smaller contributions. Interestingly, unlike the situation en¬ 
countered in semiconducting quantum wells we hnd that 
Rashba spin-orbit coupling does not lead to backscattering, at 
least to leading order. 

We have also investigated the temperature and magnetic 
held tilt angle dependence of a single-edge conductance as¬ 
suming that the main scattering sources are the scalar potential 
and the (intrinsic) spin-orbit coupling induced by the scatterer. 
In addition, the other scattering sources discussed in Sec. m 
may be generated via virtual transitions to different Landau 
levels, which makes estimating their magnitude rather diffi¬ 
cult. However, we expect them to be sub-dominant at large 
magnetic helds and therefore we have neglected them in a hrst 
approximation. 

We hope our work will shed light on the origin of the chan¬ 
nel dissipation observed in the experiment. In particular, we 
have found that disorder-induced backscattering in combina¬ 
tion with electron-electron interactions leads to a power-law 
temperature dependence of the conductance which may be ob¬ 
served experimentally. In this regard, we note that, even in the 
ferromagnet, for which the backscattering effects are milder, 
the presence of impurities, combined with repulsive interac¬ 
tions can lead to an insulating-like behaviour of the edge con¬ 
ductance. 

The model studied here can be fairly easily generalized 
to disorder potentials with more complicated spatial depen¬ 
dence. Such a study will not be pursued here. However, we 
can anticipatd^ it will yield to a different exponent in the 
temperature dependence of the conductance. A thorough un¬ 
derstanding of the origin of dissipation in the helical edges 
of graphene under a strong, magnetic field will require a 
more detailed microscopic model of the edge of the device. 
Such model should account of the existence of multiple chan- 
nelJ^and also the coupling with the gapless bulk (Goldstone) 
modes.CS 
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Appendix A: Bosonization 

In this and the following Appendices, we collect some of 
the most important results of bosonization techniqu^^ applied 
to a one-dimensional (ID) channel with two different carrier 
velocities, which is relevant to the situation studied in the main 
text. 

The Hamiltonian of the ID edge can be expressed in 
terms of the Hamiltonian density, "Hedge (a^), as Hedge = 
/ dx Hedgeix), and Hedgeix) can be fully expressed in terms 
of the fermion chiral density operators, 

Pc{x) =: tpl{x)ilJc{x) : (Al) 

where c = R,L and : ... : stands for normal order. The chiral 
densities obey the Kac-Moody algebra® 

[pc{x),pc'{x')]=^^^^^doo6{x-x'). (A2) 

Ztti 

where sn = —sl = +1- Thus, the Hamiltonian density can 
be written as® 

Hedge = 71" : [vrPr + vrPl + ‘^.QrlPrPl] :, (A3) 

where vr (vr) is the Fermi velocity of the right (left) chan¬ 
nel and guL is an effective short-range electron-electron in¬ 
teraction. It is also useful to introduce the bosonic fields 
4>c=r,l{x), which are defined by 


The new chiral fields obey a modified Kac-Moody algebra; 

[Pa{x), ppix')] = dx5{x - x ), (All) 

where a,P = ±. and = ±1. Using these commutation 
relations, the chiral densities are found to obey the following 
equations of motion: 


{dt±v±d^) p±{x,t) =Q. (A12) 


In the stationary state (i.e. for dtp± = 0), when coupled to 
two different (chiral) chemical potentials, the chiral densities 
must minimize Hedgeix) — p+p+{x) — p-p-{x), which is 
possible (upon completing the square) if we choose: 

P± = {p±{x)) = (A13) 

2Trv± 


Hence the stationary current is given by / = e(u+p+ — 
V-P-) = eK{p^ — p-)/ (27r). Note that the chemical poten¬ 
tials, p± coupling to the densities p± are physically different 
from the chemical potentials, Pr/r, coupling to the electrons. 
The latter couple to Nj^/r = f dxpfi/R{x) and are equal to 
the chemical potentials of the source and drain lead reservoirs, 
eVsjDi prespectively. However, it is possible to relate both 
sets of chemical potentials linearly,® so that the voltage bias 
V = Vs — Vr, = H(/r+ — p-)le. Hence, the quantized con¬ 
ductance for a clean ID edge channel can be recovered®® 




(A14) 


Pc{x) = ^dx(t)c[x). (A4) 

The Hamiltonian density can be brought to diagonal form: 
Hedge(a;) = ^ : [v+p\[x) + v-pi{x)\ : , (A5) 

by means of the following linear transformation; 


which is the value of the quantized conductance for a single 
channel in li = 1 units. The quantized conductance of a per¬ 
fect non-chiral Tomonaga-Luttinger liquid is the result of the 
finite contact resistance. 


Appendix B: Accounting for the Contact Resistance 


( P+{.x) 

V P-{^) 


where 


and 


^l(l + Kl-K\( pr{x) \ 
2\1-K 1 + k) \pRix) ) ’ 


K = 


vr + vl — 2 gRL 
vr + vr + 2 giiR 


(A6) 


(A7) 


In this Appendix, we shall follow the approach of Pham and 
coworkers^to account for the contact resistence of the ID in¬ 
teracting edge channel in the presence of a scatterer located at 
a: = 0. As it has been discussed in the main text, the Hamil¬ 
tonian for such system is H = Hedge + Vm = J dxH, where 
the Hamiltonian density is given by; 


H = Hedge + pRm 'iljR{x)'lpR{x) +h.c. S{x) (Bl) 


.(vr-Vr) {vr+VrY 2 

»± = ±—5— + \j[-^—) -s™- 

Note that the transformation in Eq. ( |A6| ) has the following 
properties: 

(A9) 
(AlO) 


We are interested in computing the two terminal conductance, 
which is defined from the current passing through the device, 
/ 2 t, as follows: 


_ 

~ dV 


v=o 


(B2) 


where V = Vr — Vr, is the voltage bias between the source 
and drain leads. The current l 2 t is given by the following 


p+{x) + p-{x) = pr{x) + pr{x), 
p+{x)- p-{x) = K[pRi{x) - pr{x)]. 














11 


equations (essentially Ohm’s law)® 

Vs — — = Rsht, 
e 

- - eVn = Roht, 

e 

where 

-< Ai+ + M- 
^-■ 


(B3) 

(B4) 


(B5) 


are the mean chemical potentials of the edge to the right ((2^) 
and left (/2^) of the impurity, respectively; Rs, Rd Vd) 
are the the source and drain contact resistances (voltages), re¬ 
spectively. 

In the presence of the impurity, the equations of motion for 
p± (cf Eq. |A12| l are modified and can be written as follows: 

{dt±v±dx) p±{x,t) = J±{t)5{x), (B6) 


where the operator J± = i J dx [Vm, P±[x)\ = i [Vm, iV±], 
with N± = f dx p±(x). However, note that, away from 
the scatterer (i.e. x ^ 0), the equations of motion reduce 
to Eq. ( A12| i. Thus, in the steady state, the chemical poten¬ 
tials /ij can be related to the chiral densities by means of 
Eq. ( |A13| ), that is. 




(B7) 


where are the chiral densities to the right and left of the 
scatterer. Adding equations (|B3[) and (|B4[) yields: 


lit. = 


Rq 


Vs — Vp _ 

Rs + Rd Rs + Rd 


Ib, 


(B8) 


where we have introduced Rq = Gq^ = 27r/e^ and 


Ib = -^[v+{p%- p+)+V-{p>-p<)] (B9) 


which, as shown below, is the back-scattered current. Note 
that our definition of Is differs from the definition in Ref.l40l 
by a factor of Ib can be obtained upon integraging over 
an infinitesimal interval around x = 0 the above equations 
of motion Eq. ( |B6| l and taking the expectation value, which 
yields: 

v±{pI-p<)=±{J±). (BIO) 

Hence, 

- J-)) = - ^([FM,iV+ - N_]) 

7 P 

= Nr - Nl\) 

7 g 

= --{[H,Nr-Nl]). (BID 

In deriving the last expression, we have integrated Eq. ( |A10| i 
over X and used that [Hedge, Nr/r] = 0. Next we use Nr — 


Nl = 2Nr - N and [H, A^] = 0, where N = Nr + Nr, 
which leads to 


lB = {-ie)[H,NR] = -e(^^) (B12) 

Hence, assuming a symmetric situation where Rs = Rd = 
7r/e^, it follows that 


l2t = GoV-lB, (B13) 

where Go = 12-11 is the clean channel conductance. After 

derivation with respect to the voltage bias, Eq. ( |B13| l for V = 
0 becomes: 


G = Go - SG, (B14) 

SG = (dlB/dV)v=o being the conductance across the impu¬ 
rity. The latter is computed in the following Appendix. 


Appendix C: Linear Response for the Conductance 

In this appendix, we review the calculation of the linear 
conductance of a ID channel in the presence of an impurity, 
SG, using the Kubo formula. The result is perturbative in the 
strength of the impurity potential and requires SG Go, 
which, as discussed in the main text will break down as the ab¬ 
solute temperature T —>^ 0 for repulsive interactions (K < !)■ 
We begin by introducing the backscattering current operator: 

4 = -e^ = (-ie)lH,NR] = (-ze) [A - At] , (Cl) 

where H is given by Eq. and the operator A = 

gRM4’R{0)'4’L{0)- Within linear response theory, the steady 
state current in response to a voltage bias V is (see e.g. 
Ref [36l page 561 and ff): 

/ + 00 

dt G^^t (t) (C2) 

-OO 

where the correlation function, (t) = 

—10(f)([A(f), At(0)]) where 9{t) is the Heaviside step 
function. In order to evaluate the above correlation function, 
we rely on bosonization, which uses the following representa¬ 
tion of the Eermi fields in terms of the bosonic fields (j)R/R{x) 
introduced in Appendix [A| 


= A '’**’ 

(C3) 


(C4) 


where Fc = (c = R, L) are the Klein factors, {Fc, Fc'} = 
Sec', required to ensure the anti-commutativity of the right and 
left moving Eermi fields. 

In the presence of interactions, the bosonic fields (j)RlR{x) 
do not describe the eigenmodes of the clean interacting ID 
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edge. However, we can introduce a pair of bosonic fields 
(j)± (x) reated to the chiral densities p± that diagonalize Hedge- 


P±{x) = 


(C5) 


The finite temperature correlation function of the new bosonic 
fields (for |f | ^ Ib/v±) at a; = 0 reads® 


(</'±(0,f)<^±(0,0)) = - In 


f2v±P\ . fiirt 
_ 


[ ttIb ) 


■ (C6) 


Using Eq. ( |C6| l, the correlator (f) can be evaluated and 
thus we arrive at the following expression for the current; 

I = 2e{gRMfCRe f -/(f)}, (C7) 

J —OO 


where 

f{t) = 

and the prefactor 

C = 

After integration, we obtain: 


ACe{gRM)^l3 . , ^ , teU/3 

— sin(7rAr) Im B AT + 


• (-4- ^ , Ib 

sin — ir H-sin — ir H- 

/3 \ v+J P \ v_ 


K 


-K 


( ' V 

\ 1 

\2ttIbJ 

P'^V+V- _ 


(C8) 


(C9) 


2^ 


27r 


1-2A: 


(CIO) 


where B{n,m) is the beta function. From this expression, 
after derivation with respect to U, at U = 0, Eq. follows. 
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